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INTRODUCTION 

A  semi  classical  study  presented  here  indicates  that  a  sufficiently 
short  and  intense  pulse  can  be  much  more  effective  in  inducing  a  collisional 
radiative  transition  than  cw  radiation  or  a  long  pulse,  although  the 
intensity  must  not  be  too  high  because  the  Rabl  oscillation  can  bring  down 
the  probability.  Fo»r  the  situation  of  a  molecule  physlsorbed  on  a 
crystalline  surface  and  irradiated  by  a  laser,  a  master  equation  approach, 
used  to  describe  the  time  evolution  of  the  population  of  the  vibrational 
adbond  levels,  shows  that  for  high  intensities  laser- induced  vibrational 
excitation  is  the  same  for  pulsed  and  cw  lasers. 

GAS -PHASE  MOLECULAR  COLLISIONS 

Most  theoretical  and  experimental  work  in  the  past  on  the  subject  of 
molecular  collisions  occurring  in  the  presence  of  a  laser  field  has  been 
concerned  with  the  case  in  which  the  field  consists  of  continuous  wave  (cw) 
radiation.  In  this  case  the  colliding  molecules  are  influenced  by  the  laser 
radiation  during  both  the  incoming  and  outgoing  trajectories.  If,  however, 
the  laser  field  consists  of  a  pulse  whose  duration  is  sh^jter  than  the 
collision  time  (for  a  typical  collision  time  of  10  -  10  s,  one  must 
have  a  (sub)picosecond  pulse),  only  a  part  of  the  collision  process  occurs 
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in  the  presence  of  the  laser  field.  The  outcome  of  the  collision  therefore 
may  be  significantly  different  depending  whether  the  laser  radiation  is  cw 
or  pulsed.  The  difference  shows  up  strongly  in  particular  for  collision 
processes  involving  an  electronic  tran^i£ion  that  occurs  via  curve  crossing, 
as  we  have  shown  in  our  earlier  work.  *  Below  we  give  a  brief  description 
of  such  a  transition  taking  place  in  the  presence  of  laser  radiation,  with 
special  attention  to  the  effect  of  the  temporal  duration  of  the  laser 
radiation^  An  experimental  study  of  this  effect  has  recently  been 
reported. 


Let  us  consider  an  electric-dipole  transition  induced  by  collision  that 
occurs  in  the  presence  of  laser  radiation.  Within  the  semiclassical  two- 
state  formalism,  the  Schrodinger  equation  yields  the  following  coupled 
equations  for  the  probability  amplitudes  c^(t)  and  c2(t)  for  the  two  states 
1  and  2  being  considered: 

dc. (t) 

l*:  dt-*  "llCl  +“l2c2  <la> 

dc-(t) 

dt  "  H21c1  +  H22c2  * 


where  H 


11 


and  1 
transit j 


and  H_»  represent  adiabatic  potential  curves  in  the  dressed- 
molecule**representation  for  the  two  states  1  and  2  which,  we  assume,  cross 
at  a  certain  internuclear  distance  R  ■  R  (i.e.^  H22^Rc^*Hll^Rc^  *  the 
interaction  term  H.2  is  given  by  H._  *  -  p-Ep/2,  p  is  the  transition  moment, 
is  the  electric  field  amplitude  of  the  laser  radiation.  The 
!ion  probability  1  ■»  2  is  obtained  by  solving  the  coupled  equations 
(la)  and  (lb)  with  the  initial  conditions  |c  ( 0 ) |  ■  1  and  c_(0)  ■  0.  These 
equations,  however,  cannot  be  solved  analytically  in  general  because  the 
matrix  elements  H,  vary  with  time  in  a  complicated  manner  for  a  realistic 
collision  system.  ^One  therefore  must  rely  on  either  approximate  or 
numerical  solutions.  Here  we  present  results  of  approximate  calculations 
based  on  the  Landau-Zener  method  along  with  results  of  numerical 
calculations. 


Perhaps  the  simplest  approximate  solution  to  Eqs.  (1)  is  provided  by 
the  Landau-Zener  model.  The  model  assumes  that  the  transition  is  localized 
in  a  narrow  neighborhood  of  the  crossing  point  R  over  which  H-j  *  H,, 
varies  linearly  with  R  and  H._  remains  constant.  The  transition  probability 
as  the  system  passes  through1 the  crossing  point  Rc  is  then  given  by 


P  ■  1  -  exp(-p)  (2a) 

where 

p  -  2ir[H12(Rc)]2/*vy[l  -  (|-)2]4  ,  (2b) 

c 

b  is  the  impact  parameter,  y  is  the  slope  of  I ^22*^1 1 ^  V8#  R  at  Rc'  an<*  V  *s 
the  relative  velocity  of  the  colliding  partners  at  R£.  If  the  collision 
occurs  in  the  presence  of  a  short  pulse  so  that  the  system  is  illuminated  on 
the  way  in  or  out  only,  the  transition  probability  is  given  by  P  ■  P.  On 
the  other  hand,  if  the  system  is  illuminated  both  on  the  way  in  and  out  as 
is  the  case  for  cw  radiation  or  a  long  pulse,  the  transition  probability  is 
given  by  P  -  2P(1-P).  To  exhibit  clearly  the  intensity  dependence  of  the 
transition  probability,  we  note  that  p  is  proportional  to  I  and  let  p  ■  pi. 


Then  the  short-pulse  probability  and  cw  probability  are  given  respectively 
by 


Ps  -  1  -  exp(-3I)  (3a) 

Pcw  «  2  exp(-f$I)(  l-exp(-8I)]  (3b) 

where,  for  the  case  of  a  short  pulse,  a  square  pulse  of  intensity  I  is 
assumed. 

In  Table  1  we  give  the  Landau-Zener  probabilities  P8  .and  P  for 
different  values  of  I  and  b  for  the  case  R  ■  9.5,  v  ■  1.83*10’  ,  y  «C1J.002 
and  y  »  3  (all  quantities  are  expressed  inatomic  units),  assuming  that  the 
laser  beam  is  linearly  polarized  in  the  direction  of  the  relative  velocity  v 
of  the  colliding  partners  and  that  the  transition  moment  y  is  parallel  to 
the  internuclear  axis  as  in  a  I  -  l  transition.  The  probability  P  as  a 
function  of  the  intensity  I  exhibits  the  well-known  single-peak  structure  of 
the  Landau-Zener  model,  whereas  P  monotonically  increases  with  I.  As  a 
result,  P  and  P  differ  significantly  at  high  I  for  which  81  »  1.  Thus, 
according  c£o  the8  Landau-Zener  model,  a  short  pulse  can  be  much  more 
effective  than  cw  radiation  or  a  long  pulse  in  inducing  a  transition,  if  the 
pulse  intensity  is  sufficiently  high. 

Table  1.  The  cw  probability  P  and  short-pulse  probability  P  calculated 
using  the  Landau -Zener  model.  8 


b(a.u. ) 

I(W/cm2)  1q8 

5*10® 

109 

5*109 

io10 

5*1010 

10U 

5*10U 

P 

0.053 

0.23 

0.37 

0.39 

0.13 

0.00 

0.00 

0.00 

1 

cw 

P 

s 

0.027 

0.13 

0.24 

0.74 

0.93 

1.00 

1.00 

1.00 

p 

0.045 

0.20 

0.33 

0.43 

0.18 

0.00 

0.00 

0.00 

5 

cw 

p 

0.023 

0.11 

0.21 

0.69 

0.90 

1.00 

1.00 

1.00 

s 

In  order  to  check  the  accuracy  of  the  data  presented  in  Table  1,  we 
also  have  calculated  the  transition  probability  by  numerically  integrating 
Eqs.  (1).  For  this  calculation,  we  have  assumed  that  the  potential  curves 
satisfy 

H22  ’  H11  "  3  exp(-0.737R)  -  0.0027  ,  (4) 

which  yields  R  ■  9.5  and  y  ■  0.002.  We  have  also  assumed  a  straight-line 
constant  velocity  ("^  v)  trajectory  for  simplicity  of  calculation  and  have 
chosen  v  ■  1.83  *  10  and  y  ■  3,  as  before.  The  result  of  our  calculation 
is  summarized  in  Table  2.  Comparison  of  Tables  1  and  2  indicates  that  the 
Landau-Zener  model  yields  relatively  accurate  probabilities,  at  low 
intensities,  but  the  model  fails  at  high  intensities  (I  >  lO11  W/cni  ).  This 
is  mainly  due  to  the  inability  of  the  Landau-Zener  model  to  correctly 


describe  the  Rabi  oscillation  of  the  molecular  system  at  a  high  laser 
intensity.  Since  Ps . .  decreases  due  to  the  Rabi  oscillation  as  I  is 
increased  beyond  -10  W/cm  ,  there  is  a  finite  range  of  the  intensity  over 
which  a  short  pulse  is  much  more  ef fective-than  cw  radiation.  For  the  model 
system  being  considered,  this  range  is  10IU  -  10U  W/cm. 

Table  2.  The  cw  probability  P  and  short-pulse  probability  P  calculated 
by  numerical  integratiSn  of  Eqs.  (1).  s 


P  0.029  0. 14  0.27  0.73  0.88  0.92 

s 


0.034  0.17  0.28  0.41  0.13  0.13 

2  cw 

P  0.019  0.091  0.17  0.71  0.93  0.93 

s 


Ps  0.86  0.73  0.68  0.60  0.61  0.54 


P  0.23  0.40  0.46  0.49  0.48  0.47 

5  cw 

Ps  0.87  0.72  0.64  0.58  0.61  0.63 


In  view  of  the  recent  experimental  study  of  the  short-pulse  effect 
considered  here,  we  show  in  Table  3  the  total  number  N  of  transitions  per 
target  atom  per  unit  pulse  energy,  an  experimentally  observable  quantity, 
for  different  values  of  the  pulse  duration  and  Intensity.  Although  this 
quantity  was  obtained  by  integrating  the  Landau-Zener  probability  over 
impact  parameter,  the  data  shown  in  the  table  is  expected  to  be  accurate 
since  the  Landau-Zener  model  yields  accurate  probabilities  at  the  values  of 
I  considered  here.  As  can  be  seen  from  the  table,  thera^is  very  little 
difference  between  long  and  short. -pulses  at  10°  -  10y  W/cm.  The 
difference,  however,  shows  up  at  10*  W/cm  ,  as  a  noticable  jump  in  is 
seen  in  going  from  10  to  1  ps  (the  collision  time  of  the  system  being 
considered  is  -  2.5  ps). 


Table  3. 

Total  number  of  transitions 
(in  arbitrary  units). 

per  target 

atom 

per  unit  pulse  energy 

I(W/cm2) 

\  pulse 

\  duration  (ps) 

1000 

100 

10 

1 

0.1 

108 

19.8 

19.8 

19.9 

20.2 

20.2 

109 

15.1 

15.1 

15.7 

18.5 

18.5 

1010 

1.45 

1.57 

2.81 

8.75 

8.84 

In  conclusion,  the  results  of  our  calculations  summarized  in  Tables  1, 
2  and  3  indicate  that  a  short  pulse  can  be  much  more  effective  in  inducing  a 
collisional  radiative  transition  than  cv  radiation  or  a  long  pulse.  This 
short-pulse  effect  becomes  significant  if  the  pulse  intensity  is 

sufficiently  high  that  the  probability  of  transition  at  the  crossing  point 

is  high  (-  1).  On  the  other  hand,  the  intensity  must  not  be  too  high 

because  the  Rabi  oscillation  can  bring  down  the  probability  to  a  low  value 

at  an  extremely  high  intensity. 

ADBOND  DYNAMICS 

Level  Populations 

We  now  turn  to  consider  a  molecule  physisorbed  on  a  crystalline  surface 
and  irradiated  by  a  laser.  Due  to  the  weak  van  der  Waals  bond,  the 
admolecule  has  a  series  of  vibrational  states  for  the  motion  perpendicular 
to  the  surface.  Under  certain  conditions  it  is  reasonable  to  assume  that  a 
laser,  with  frequency  w.  ,  is  in  reasonance  with  only  one  pair  of  these 
levels,  say  |g>  and  |e>.  The  transition  frequency  is  u  ■  w  -  u  >0;  the 
detuning  is  defined  as  A  »  u.  -  w  .  Relaxation  of  theefdbonS  occSrs  due  to 
its  interaction  with  the  lattice  vibrations  of  the  substrate.  In  the  Markov 
approximation,  this  is  described  through  the  occurrence  of  rate  constants 
a  ,  in  a  set  of  first-order  differential  equations  for  the  level  populations 
ana  coherences  of, the  reduced  density  operator  of  the  adbond,  the  so-called 
master  equation. 

For  pulsed  lasers  with  a  pulse  duration  At  short  compared  to  the 
relaxation  times  (■  the  inverse  rate  constants)  of  the  adbond,  the  latter 
can  be  neglected  during  the  pulse.  Then,^he  time  evolution  during  the 
pulse  is  given  by  the  optical  Bloch  equations, 

-  AR2(t) 

-  -AR.(t)  +  n  (t)R^(t)  (5) 

1  P  J 

-  -0  (t)R-(t)  , 

P  2 


dRjU) 

dt 

dR2(t) 

dt 

dR3(t) 

dt 


1 


\ 


V\ 


where  fl  (t)  ■  y.fg(t)/h  is  the  time-dependent  Rabi  frequency;  2(t)  ■ 
£_{t)cos\wLt)  is  £ne  electri^  field  amplitude  of  the  laser,  with  slowly- 
varying  envelope  iL(t);  and  y  again  is  the  transition  dipole  moment.  The 
real-valued  quantities  R.  are  defined  by  the  populations  of  and  coherences 
between  the.Qtwo  adbond  levels  coupled  by  the  laser,  in  the  rotating  frame, 
according  to1 


R 

R 

R 


1 

2 

3 


P  +  P 
8®  - 

*i  (P„ 

ge 

P  -  P 
e  g 


eg  - 
-  P 


(6) 


s 


The  total  exciting  effect  up  to  time  t  of  a  laser  is  characterized  by 


e(t)  -  /c  dt'  n  ( t * ) 

ta>  p 


(7) 


Here,  we  are  only  interested  in  the  effect  of  a  complete 
that  the  integral  can  be  taken  over  a  single  pulse 
independent  parameter  0 .  Let  R° ,  R°  and  R^  be  the  values 

pu^ge.  Then,  for  zero-detuning,  the  values  right  after  the 
byIU 


laser  pulse,  so 
to  give  a  time- 
just  before  the 
pulse  are  given 


R^At)  -  RJ 

R2(At)  -  R“cos(0)  +  R°sin(8)  (8) 

R^CAt)  «  R“cos(0 )  -  R“sin(8)  , 


where  the  start  of  the  pulse  is  taken  as  zero  point  of  time,  and  At  is  the 
pulse  duration.  We  shall  only  consider  the  situation  of  maximal  excitation, 
which  is  obtained  for  8  -  it,  the  so-called  n-pulses. 

Instead  of  a  single  pulse,  we  shall  consider  a  series  of  equally  spaced 
ir-pulses,  with  interval  time  t  and  pulse  duration  At  <<  t  .  Assume  that 
the  system  is  initially  in  thermSl  equilibrium,  i.e.,  R^  »  0?  R2  ■  0  and  R^ 
•  R-j(eq).  Then,  after  a  ir-pulse,  and  R-  are  zero  again,  but  R,  is 
changed  into  its  opposite:  R^(At)  ■  -  R^Ceq;.  Between  two  consecutive 
pulses  the  adbond  evolves  in  time  through  its  relaxation  against  |hg  lattice 
vibrations.  This  process  can  be  described  by  a  master  equation.  In  the 
present  situation  only  the  populations  have  to  be  considered.  However, 
unlike  the  coupling  with  the  laser,  transitions  to  and  from  all  vibrational 
states  are  possible.  With  the  result  of  the  laser  pulse  as  initial 
condition,  the  formal  solution  of  the  master  equation  is  (At  <  t  <  t^), 

P(t)  -  P(eq)  +  e  Wt{P(At)  -  P(eq)}  ,  (9) 

**  *m  am 


where  W  is  the  matrix  formed  by  the  transition  rate  constants  a^  and  has 
elements 


£(t)  is  the  vector  formed  by  the  populations  P  (t),  end  P(eq)  refers  to  the 
thermal  equilibrium  distribution,  given  by  W  P(eq)  ■  0. 

After  a  number  of  pulses  have  passed,  the  adbond  will  reach  a  quasi¬ 
steady  state,  wherein  the  time  evolution  of  the  populations  is  the  same  in 
each  interval  t  (see  Fig.  1).  Then,  combining  Eqs.  (8)  and  (9),  two  sets 


Fig.  1.  Schematic  drawing  of  the  time  evolution  of  the  level  occupations  for 

a  pulse  sequence  in  the  quasi-steady  state.  Curve  a  represents  R^; 

curve  b  gives  P  +  P  ;  and  curve  c  gives  P  for  any  other  level  n  * 
8  e  n 

e  or  g.  * 

of  equations  are  obtained,  which  together  determine  the  ^nstants  Pn(At)  in 

terms  of  the  rate  constants  a  .  and  the  interval  time  t  : 

ntc  p 


-Wt 

P(t  )  -  P(eq)  +  e  p  (P(At)  -  P(eq)} 

P  (t  )  -  P  (At)  -  P»  (n  *  e,g) 
n  p  n  n 

P  (t  )  -  P  (At)  -  P» 

g  p  e  g 

p_(t_)  -  p„ut)  -  p;  . 
e  p  g  e 


(11) 


(12) 


For  a  continuous  wave  laser,  the  Rabi- frequency  is  a  constant,  0  . 
Now,  the  relaxation  processes  cannot  be  neglected  during  the  interaction 
with  the  laser.  The  rate  constants  e^  give  rise  to  additional  terms  in  the 


a 


optical  Blocb.equations ,  in  a  way  analogous  to  their  occurence  in  the  master 
equation.  ’  After  some  time  a  true  steady  state  will  be  reached,  in 
which  the  populations  are  determined  by 

l  *«kY->  •  I  vV-*  <n  *  ••«> 

k  k 

I VY*' '  l  Vk<-> '  2  o "h?  <Y*>  •  Y*»  (13> 

k  k  2*  eg+^a 

l  *,kp.<">  -  I  *k.V->  +  5  1^-2  <%<-'  ‘  P.<-»  • 

k  k  2  eg+za 

where  *  ^®ek  +  *gk^'  an<*  t*ia  steady-state  values  are  denoted  by  Pn(®). 


Pulsed  Laser  Versus  Continuous  Wave  Laser 

To  compare  the  effect  of  a  pulse  train  with  the  effect  of  a  continuous 
wave  (cw)  laser  on  the  dynamics  of  some  system,  a  criterion  is  needed  to 
compare  both  lasers.  Below  two  possible  criteria  will  be  defined.  First, 
it  is  required  that  both  lasers  have  equal  average  power.  Because  the  Rabi- 
frequency,  n(t),  is  proportional  to  the  electric  field  of  the  laser,  the 
intensity  is  proportional  to  CT(t).  Performing  the  time-average  for  a 
series  of  u-  ulses  then  gives: 

°cw  •  *2/<it  y  •  (i4> 

A  second  possible  criterion  is  to. require  the  average  dissipation  from  both 
lasers  to  be  equal.  This  leads  to1 

<P.  .  pM/t  -  (n?  /r  HP  (-)  -  P  (-))  ,  (15) 

g  e  p  cw  eg  g  e 

where  P°,  P#,  P  (<*>)  and  P  (•)  follow  from  Eqs.  (11)-(13). 

g  e  g  e 

The  efficiency  of  a  pulsed  laser  over  a  cw  laser  to  induce  some  process 
(e.g.,  desorption  or  chemical  reaction)  will  be  different  for  each  process. 
Here  we  restrict  ourselves  to  the  effect  on  the  population  of  '"he  upper  one 
of  the  two  levels  coupled  by  the  laser  for  a  two-level  system.  A  more 
extensive  discussion  is  given  elsewhore.  Let  us  define  the  efficiency  c 
as  the  ratio  of  the  average  excited  level  population  P  (av)  under  pulsed 
laser  irradiation  to  the  steady  state  population  P  ?»)  under  cw  las^y 
Irradiation.  For  a  two-level  system,  Eqs.  (11)-(13)  are  easily  solved. 

We  are  interested  in  the  situation  where  the  laser- induced  population  is 
substantially  larger  than  the  equilibrium  thermal  occupation  of  the  levels. 
Then,  neglecting  the  latter,  the  efficiency  is1 


(1  +  e  rtp) 


[1  -  e 

rt 


1  2  2 

+  n 

2  cw 


For  a  given  integrated  intensity  ot  one  pulse,  the  average  intensity  of  the 
pulse  train  is  proportional  to  t^  .  Then  a  definition  for  a  strong/weak 


pulsed  laser  is  £tp  $  l *  Analogously  a  strong/weak  continuous  wave  laser 
is  defined  by  (Gcw/r  )  being  much  larger/smaller  than  one. 


Using  these  definitions,  it  follows  from  Eq.  (16)  that  c  •  1  in  the 
strong  laser  limit.  This  can  easily  be  understood  when  it  is  realised  that 
both  for  the  pulsed  laser  and  the  cw  laser  the  system  is  then  continuously 
driven  into  saturation  (R^  »  0).  Using  criterion  1  (Eq.  (14))  leads  to 


i  *rt 
1  LL  P 


2  2 

r  +  2ir 

cw 


1  + 


■Yt 


rat 

_2 


(17) 


2 

In  the  weak  laser  limit  this  gives  c  ■  FAt/x  «  1,  and  thus  c  can  be  made 
arbitrarily  small  by  decreasing  the  pulse  duration  At.  Thj  explanation  is 
that  the  intensity  of  the  pulsed  laser  is  proportional  to  G*At,  whereas  the 
exciting  effect  is  only  proportional  to  G  At.  Decreasing  A?,  while  G^At  is 
held  constant,  will  decrease  G  At  too.  ^For  a  cw  laser  both  the  intlnsity 
and  the  exciting  effect  are  propo?tional  to  fl2  .  The  latter  can  be  regarded 
as  the  limiting  form  of  a  pulsed  laser  with  $¥  -+  <*>,  and  therefore  will  give 
the  most  efficient  form  of  excitation. 


For  criterion  2,  c  *  1  is  obtained,  independent  of  the  laser  power. 

This  result  reflects  the  fact  that  the  energy  flow  into  the  substrate  (is 

absorbed  energy)  is  proportional  to  the. average  occupation  of  the  excited 

level  and  not  on  how  P  varies  in  time.  * 

e 
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